Gravitational wave energy spectrum of a parabolic encounter 
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We derive an analytic expression for the energy spectrum of gravitational waves from a parabolic 
Keplerian binary by taking the limit of the Peters and Mathews spectrum for eccentric orbits. This 
demonstrates that the location of the peak of the energy spectrum depends primarily on the orbital 
periapse rather than the eccentricity. We compare this weak-field result to strong-field calculations 
and find it is reasonably accurate (~ 10%) provided that the azimuthal and radial orbital frequencies 
do not differ by more than ~ 10%. For equatorial orbits in the Kerr spacetime, this corresponds 
to periapse radii of r p > 20M. These results can be used to model radiation bursts from compact 
objects on highly eccentric orbits about massive black holes in the local Universe, which could be 
detected by the Laser Interferometer Space Antenna (LISA). 
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I. INTRODUCTION 

An important source of gravitational waves for the 
proposed space-based gravitational wave detector, the 
Laser Interferometer Space Antenna (LISA) [l], are 
the inspirals of stellar-mass compact objects into mas- 
sive black holes in the centres of galaxies. During the 
last few years of inspiral these systems generate continu- 
ous gravitational waves in the LISA band, which will al- 
low the detection of as many as several hundred systems 
out to redshift z ~ 1 ||. However, prior to this phase, 
the inspiraling object spends many years on a highly ec- 
centric orbit, generating bursts of gravitational radiation 
at each periapse passage. LISA could resolve individual 
bursts from sources in the nearby Universe. Initial es- 
timates [|| suggested a LISA event rate of ~ 18 yr _1 , 
including ~ 15 yr _1 from the centre of the Milky Way, 
although this was subsequently revised downwards to 
~ 1 yr _1 with negligible contribution from extragalac- 
tic sources If even a single burst from the Galactic 
centre is detected during the LISA mission, this will pro- 
vide an unparalleled probe of the structure of spacetime 
there. 

The spectrum of radiation from these bursts will be 
well approximated by the spectrum of a parabolic orbit 
In this note we derive an analytic approximation to this 
spectrum by taking the limit of the Peters and Math- 
ews [(| 0] (PM) energy spectrum for eccentric Keplerian 
binaries. We show that the peak of the spectrum de- 
pends primarily on the orbital periapse and only weakly 
on the eccentricity. We also estimate the range of va- 
lidity of the approximation (in Sec. by comparing 
to numerical Teukolsky data, finding that it is a good 
approximation for equatorial orbits in Kerr with peri- 



apse r p > 20M. The parabolic spectrum takes a neat 
analytic form; deriving it from the bound spectrum will 
allow corrections for high-eccentricity bound orbits to be 
found in the future. We hope this note will be a useful 
resource for future work on gravitational radiation from 
high-eccentricity orbits. 



II. PARABOLIC LIMIT 

A. Energy Spectrum 

For an orbit of eccentricity e with periapse radius r p , 
Peters and Mathews p give the power radiated into the 
nth harmonic of the orbital angular frequency as 
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where the function g(n, e) is defined in terms of Bcsscl 
functions of the first kind 
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The Keplerian orbital frequency is 
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1 We use "parabolic" to refer to marginally bound orbits. 
Marginally bound Keplerian orbits are parabola; in curved space- 
times they do not retain such a simple shape. 



where lj c is defined as the angular frequency of a circular 
orbit of radius r p . The energy radiated per orbit into the 
nth harmonic, that is, at frequency uj n = nui, is 
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as e — !• 1 for a parabolic orbit, u>i — > as the orbital 
period becomes infinite. The energy radiated per orbit is 
then the total energy radiated. The spacing of harmonics 
is Aw = lo\ , giving energy spectrum 



dE 
did 



Wi = E(n). 



(5) 



Changing to linear frequency 27r/ = ui, 
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where the function £(n, e) is defined in the last line. For 
a parabolic orbit, we must take the limit of £(n, e) as 
e-> 1. 

We simplify ^(n, e) using the recurrence formulae 
(Watson || 2.12) 
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and eliminate n using 
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where / = uj n juj c — f n /f c is a dimensionless frequency. 
To find the limit we define two new functions 
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To give a well-defined energy spectrum, both of these 
must be finite. 

The Bessel function has an integral representation 
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we want the limit of this for v — ¥ oo, z — > oo, with z < v. 
Using the stationary phase approximation, the dominant 
contribution to the integral comes from the regime in 
which the argument of the cosine is approximately zero 
(Watson H 8.2, 8.43), for small i9: 

JJz)~- [ cos ( vd - zi) + -tf 3 ) d-d (13) 
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this last expression is an Airy integral and has a standard 
form (Watson U 6.4) 
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where A„(z) is a modified Bessel function of the second 
kind. Using this to evaluate the limit gives 
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For our case, 
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and the first limiting function is well defined, 
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To find the derivative we combine © and (fTB"]) . and 
expand to lowest order yielding 
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We may re express the derivative using the recurrence 
formula (Watson [8] 3.71) 
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And so finally we obtain the well-defined 
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Having obtained expressions for A(f) and B(f) in 
terms of standard functions, we can now calculate the 
energy spectrum for a parabolic orbit. From ([7]) 
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where we have used the limit 
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FIG. 1. The relative energy (per orbit) spectrum l(n,e) for 
e = 0.2 (heavy line), e = 0.5 (medium line), e = 0.7 (light 
line), and the limiting result for e = 1 (dashed line) versus 
frequency. 



This agrees with the e = 1 form of Turner's result, which 
was computed by direct integration along unbound or- 
bits [1]. Figure Q] shows how i(n,e) changes with ec- 
centricity including our result for a parabolic encounter 
(cf. Figure 3 of Peters and Mathews @). Although 
more power is radiated into higher harmonics, the peak 
of the spectrum does not move much: it is always be- 
tween f = f c and / = 2/ c , with / = 2/ c for e = and 
/ ~ 1.637/ c for e = 1. 



B. Total Energy 

To check the validity of this limit we can calculate the 
total energy radiated by integrating (|23l) over all frequen- 
cies, or by summing the energy radiated into each har- 
monic. These must yield the same result. Summing: 
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where we have used equations ([T]), (j3|) and (|4|. Peters 
and Mathews [6j provide the result 
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Using this, 
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which is perfectly well behaved as e — > 1, 
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Integrating the energy spectrum (l23|) gives 
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The integral can be evaluated numerically as 

r°° a o ^ 

J o 1(f) df = 12.5216858...= 

The two total energies are consistent, E- mt = E s 
III. APPLICABILITY 



A. Limit of Approximation 

The PM approach assumes Keplerian orbits in flat 
spacetime. This should be a valid approximation in 
the weak-field regime far from a massive body. To find 
the limit of this approximation, we can compare the 
PM results with those from more accurate techniques. 
Energy spectra for parabolic orbits do not seem to be 
available in the literature yet, so we will make do with 
the total energy fluxes calculated by Martel [10(, who 
uses time-domain black hole perturbation theory for a 
Schwarzschild black hole of mass M. Figure [2] shows the 
ratio of the two energies as a function of periapse dis- 
tance. As expected the PM result is more accurate for 
larger periapses. The agreement worsens as the periapsis 
decreases. At r p = AM, corresponding to the radius of 
the innermost stable circular orbit (ISCO), the energy 
flux calculated by Martel diverges, so the ratio tends to 
zero. This divergence is because in Schwarzschild (or 
Kerr) spacetime a parabolic orbit may have a zoom- whirl 
structure where it undergoes a number of near circular 
rotations (whirls) about the black hole. As the radius of 
the ISCO is approached, the number of whirls tends to 
infinity (in the absence of radiation reaction), so an in- 
finite amount of energy is radiated. Figure [2] also shows 
how the ratio of energies follows the number of rotations, 
defined as N = A</>/27r, where A</> is the total change in 
the azimuthal angle over one orbit. As N increases, the 
PM approximation worsens because the Keplerian orbit 
does not include this extra rotation. The accuracy of 
the PM result deteriorates rapidly once the orbit transi- 
tions to a zoom-whirl trajectory and is therefore far from 
parabolic in shape. 

The PM result is accurate to ~ 10% for orbits with 
N < 1.1. We will adopt this as a cutoff point. For an 
equatorial orbit in Kerr spacetime, 
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FIG. 2. Ratio of the total energy radiated as calculated using 
the Peters and Mathews [y] approach to that calculated by 
Martel [l(| using black hole perturbation theory (solid line) 
versus periapse radius r p . The latter approach should give 
more accurate results. Also shown is the reciprocal of the 
number of rotations 1/N (dashed line). The Keplerian limit 
corresponds to N = 1. 



L z is the specific angular momentum about the z-axis; a 
is the spin parameter, and we have adopted units with 
G = c = 1. We will find it useful to define 



r± = M ± ^M 2 - a 2 , 



and the two nonzero roots of the cubic 
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the periapsis is the larger root r p > r\. This equation 
implicitly gives L z as a function of r p . The integral may 
be rewritten as 
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where 
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This may be evaluated using elliptic integrals (Grad- 
shteyn and Ryzhik [llj] 3.131.8, 3.137.8) 
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where II(n|m) = J* di?/(l — rtsin 2 d) y/l — m sin 2 i9 is 
the complete elliptic integral of the third kind. In the 
limit of a — > we recover the Schwarzschild result [13 
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FIG. 3. Periapse radius corresponding to N = 1.1 as a func- 
tion of spin parameter a (solid line). The curve may be ap- 
proximated by a straight line r v — —3.91a + 17.36M (dashed 
line) . 



where K{m) = f* /2 d$/Vl — m sin ■& is the complete 
elliptic integral of the first kind. Figure [3] shows the peri- 
apsis for which N = 1.1 for a range of spins. Equatorial 
orbits with larger periapses should be reasonably approx- 
imated by the PM result. 

Nonequatorial orbits are more complicated because of 
the additional precession of the orbital plane. This extra 
rotation will mean that the PM approach is less accurate; 
however, this should be subdominant to the perihelion 
precession effect and so the cutoff periapsis should not 
be much larger than for the equatorial case. 



B. Astrophysical Implications 

Considering bursts from the Galactic centre, orbits 
with periapses of r p < 120M could generate bursts that 
would be detectable with LISA 0, It is therefore 
likely that any such burst that was detected would be 
in the regime of validity of the Peters and Mathews ap- 
proach, r p > 20M for equatorial orbits. The results de- 
scribed in this note will therefore have application in that 
context, and it should be possible to explore the major- 
ity of parameter space using this approximation. The 
most interesting orbits, those which come deep within 
the strong-field region of the black hole's spacetime, will 
be beyond the range of validity of this approximation, 
but these represent a small subset of all plausible events. 

This result may also be applicable for studying 
parabolic encounters between stellar mass black holes; 
these may occur in densely populated environments such 
as globular clusters [l]| or the Galactic centre [Tl j. 
Bursts from these encounters should be detectable with 
near-future ground-based detectors, such as the Ad- 
vanced Laser Interferometric Gravitational- Wave Obser- 
vatory [HEii. 
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